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ABSTRACT 

Context. Accurate star formation histories (SFHs) of galaxies are fundamental for understanding the build-up of their stellar content. 
However, the most accurate SFHs - those obtained from colour-magnitude diagrams (CMDs) of resolved stars reaching the oldest 
main sequence turnoffs (oMSTO) - are presently limited to a few systems in the Local Group. It is therefore crucial to determine 
the reliability and range of applicability of SFHs derived from integrated light spectroscopy, as this affects our understanding of 
unresolved galaxies from low to high redshift. 

Aims. To evaluate the reliability of current full spectral fitting techniques in deriving SFHs from integrated light spectroscopy by 
comparing SFHs from integrated spectra to those obtained from deep CMDs of resolved stars. 

Methods. We have obtained a high signal-to-noise (S/N ~ 36.3 per A) integrated spectrum of a field in the bar of the Large Magellanic 
Cloud (LMC) using EFOSC2 at the 3.6 meter telescope at La Silla Observatory. For this same field, resolved stellar data reaching the 
oMSTO are available. We have compared the star formation rate (SFR) as a function of time and the age-metallicity relation (AMR) 
obtained from the integrated spectrum using STECKMAP, and the CMD using the lAC-star/MinnlAC/IAC-pop set of routines. For the 
sake of completeness we also use and discuss other synthesis codes (STARLIGHT and ULySS) to derive the SFR and AMR from the 
integrated LMC spectrum. 

Results. We find very good agreement (average differences ~ 4.1 %) between the SFR(t) and the AMR obtained using STECKMAP 
on the integrated light spectrum, and the CMD analysis. STECKMAP minimizes the impact of the age-metallicity degeneracy and has 
the advantage of preferring smooth solutions to recover complex SFHs by means of a penalized We find that the use of single 
stellar populations (SSPs) to recover the stellar content, using for instance STARLIGHT or ULySS codes, hampers the reconstruction 
of the SFR(t) and AMR shapes, yielding larger discrepancies with respect to the CMD results. These discrepancies can be reduced if 
spectral templates based on known and complex SFHs instead of just single SSPs are employed. 

Key words, galaxies: stellar content — galaxies: formation — galaxies: individual (Large Magellanic Cloud) — galaxies: photometry 
— techniques: spectroscopic — methods: observational 


1. Introduction 


The study of the star formation histories (SFH) of galaxies is 
a key element in understanding their past evolution. How the 
baryonic component of a galaxy has formed and evolved should 
be reflected in its stellar populations. Due to observational con¬ 
straints, different methodologies are applied to the study of the 
stellar content in galaxies (see |Peletier|[2013[ and references 


therein). Intrinsically, these methodologies are very different and 
are affected by different sources of errors. 

Deep Colour-Magnitude Diagrams (CMD) reaching the old¬ 
est main sequence turnoff (oMSTO) are generally regarded as 
the most direct and reliable observables in order to obtain a de¬ 
tailed SFH of a galaxy ( [Gallart et al.||2005| ). This is because at 
magnitudes brightei[^than the oMSTO, stars along the main se¬ 
quence are distributed in a sequence of age: short lived, young, 


Based on observations obtained at the 3.6m ESO telescope on La 
Silla (Chile) and with the Hubble Space Telescope, operated by NASA. 


^ Below the oMSTO, there is basically no age information, but the 
luminosity function of the low main sequence can be used to obtain 
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massive stars are the bluest and brightest, and less massive and 
therefore longer lived stars are progressively fainter and redder. 
While there remains some age-metallicity degeneracy in the po¬ 
sitions of the stars on the main sequence, it is significantly less 
than that seen in other stellar evolutionary phases such as the red 
giant branch (RGB), or the horizontal-branch. In such features, 
the position of the stars are mainly determined by metallicity 
and there is very little age sensitivity. By combining the infor¬ 
mation of the position of the stars in the main sequence of a 
CMD reaching the oMSTO, with the number counts across it, 
it is possible to minimize the remaining age-metallicity degen¬ 
eracy and obtain accurate, detailed and reliable SFHs, including 


Dolphin] 

2002l|Cole et al.|2007[|Noel et al.|2009 |Monelli et al. 

2010a 

£] 

Hidalgo et al.|201 l[[de Boer et al.|2012l Meschin et al. 


so far limited to a few dozens of nearby systems, those within 
a distance of ^ 1-2 Mpc (see |McConnachie et al.||2005l for an 


updated census of the Local Group, and distance of its members 
to the Milky Way). The wide variety of morphologies and char¬ 
acteristics found in galaxies forces us to study systems at larger 
distances, where the analysis of individual stars is unfeasible. 
In these systems, because of the limited spatial resolution, we 
need to derive the stellar content using colours or spectr al infor- 
mation coming from integrate d stellar populations (e.g.,|de Jong| 


1996 


[MacArthur et al.|2QQ4[|Perez & Sanchez-Blazquez|2011t 


Sanchez-Blazquez et al.|2011| among many others) 


In order to characterise the composite stellar populations 
of unresolved galaxies, and in particular, of galaxies at high 
redshift, significant effort has been put since the early 80’s 
to interpret integrated stellar populations (e.g., |Pickles|[l985t 


Bica||1988) ). Broad band colours were first used to separate 
between old/young, metal poor/rich stellar populations dSearle 


[eTaLlfTWSl [Tinsley & Gunn|[T?7^ jPeletier & Balcells|fT996 ) 
However, this type of analysis is affected by the known age- 
metallicity degene racy and the effect of dust (e.g., j Worthey j 1994 


[Ganda et al.|2009| ). Spectroscopic analysis can help to minimize 
that degeneracy. The very first spectroscopic approach was de¬ 
veloped by|de Vaucouleurs & de Vaucoul eurs|dl959| ), later fol¬ 
lowed by [5pinrad(^1962| ) and jAlloin et al.(^197lt . The use of line 
equivalent widths or line-strength indices, taking into account 
their dependence on metallicity and age, became popular in the 


80’s and 90’s (e.g.. Rose] 19841 [Faber et al. 

19851 [Bica & Alloin 

1986a|bl Bica [19881 Gorgas et al.[ 1993[ 

Worthey||1994 

[Bica 

et al.[ 1994[ [Vazdekis et al. [19961 [Worthey & Ottaviani 

19971 

Kauffmann et al.|2003[). Indices have been used to obtain single 


information on the initial mass function (IMF) of low mass stars (e.g., 
Kalirai et al.|2013| l 


used to study the actual shape of the SFH in galaxies with com¬ 
plex stellar compositions. 

With the goal of enlarging the spectral coverage (and thereby 
maximizing the information used from the observed spectra) 
and of being able to reproduce the SFH of complex systems, 
spectral energy distribution (SED) fitting codes have been devel¬ 
oped such as MOPED ([Heavens et al.|2000[|R5chardt et al.|2001|). 


VESPA ( [Tojeiro et al. 2007[ ) or STARLIGHT ( [Cid Fernandes et a 


[2005[ ) along with other approaches. In this approach, the overall 
shape of the observed spectrum is fitted through a combination 
of stellar population models. These codes are sensitive to prob¬ 
lems in the data such as flux calibration or extinction errors, as 
they also take into account the continuum in the fit. At t he same 
time, the first so-called “full spectrum fitting codes” (e.g., [Qcvi^ 


[et al.|2006b|a|[Koleva et al.pOO^l STECKMAP and ULySS re¬ 


spectively) became available. These techniques avoid problems 
in the flux calibration and extinction in the spectra by fitting a 
polynomial to the shape of the continuum. Both SED and full 
spectrum fitting codes are better at reducing the impact of the 
age-metallicity degeneracy than line-strength indices analysis as 
they maximize the information used from the observed spectrum 
([Koleva et al. [2008] [Sanchez-Blazquez et al. [20 lT| [Beasley et al.[ 

[2015[ ). Along with this evolution in the fitting codes, there has 
been a huge progress in the modeling of stellar populations (e.g.. 


Bruzual & Charlot|2003l[Lee & Worthey |2005l[Schiavon|2007l 

Conroy et al.|2009l 

Vazdekis et al.|2010[) partly based on the im- 

provement of stellar libraries (e.g., Prugniel & Soubiran[[2001[ 

Le Borgne et al.[[2003l Valdes et al.[[2004l Sanchez-Blazquez[ 

et al.[[2006t 

Prugniel et al.[[2007[), isochrones and evolutionary 

tracks (e.g.. 

Girardi et al.[[2000al[Pietrinferni et al.[[2004|[Bres-[ 

[san et al. 

2012|[Pietrinferni et al. 

20131), and studies on the IMF 

([Salpeter 

[1955[ [Vazdekis et al.[ 

|Kroupa||20011 Weidner 


Various workers have tested the various SE D and full spec- 
trum fitting codes using artificial spectra (e.g., [Cid Fernandes 


et al.[[2005|[Qcvirk et al.[[2006a[ [Sanchez-Blazquez et al.[[2011[ 


stellar population (SSP) equivalent values (age and metallicity) 
or even to derive the whole SFH shape (e.g., [Proctor et al.[2000[ ). 
The use of indices has been relatively successful in the char¬ 
acterization of stellar populations in “simple” systems such as 
globular clusters or ellip tical galaxies (e.g., [Peletier et al.|2007[ 
[Kuntschner et al.|2010[ ). Yet, these studies i) are still hampered 
by the age-metallicity degeneracy hindering the derivation of a 
reliable SFH, ii) are biased towards the youngest stellar popula¬ 
tions which contribute much to the light while involving a small 
mass fraction, iii) make use of a limited part of the observed 
spectra, and iv) are quite limited when replicating the whole 
shape of the SFH. Combinations of different indices allow us 
to better constrain the stellar information, but this approach is 
still affected by the above limitations. Other approaches must be 


Koleva et al. 2011). In particular, Koleva et al. (2011J test 
STARLIGHT, STECKMAP, and ULySS with 48 mock spectra with 
different known SFHs. A general result of these studies is that 
the final success in the recovery of the stellar population content 
lies in the quality of the spectrum: the better the S/N and the 
resolution of the observed spectra, the more reliable the stellar 
population determination will be. 

To further test the reliability of full spectrum fitting tech¬ 
niques, it is of crucial importance to thoroughly compare the 
(in principle) more reliable CMD results to those obtained us¬ 
ing these codes with actual data. This exercise should reveal 
where they fail and where they succeed, and thus inform about 
the spectral ranges and features which are more likely to improve 
these techniques. Such analysis is not straightforward. Not only 
do CMDs containing stars as faint as those on the oMSTO need 
to be studied, but also high resolution and high S /N integrated 
spectra must be obtained. The latter is hampered by the often low 
surface brightness of nearby resolved systems. Although previ¬ 
ous studies have tried to do something similar in star clusters 


(e.g., 

Gibson et al.| 19991 [ 

Beasley et al.|2002| de Grijs & Anders 

2006 

Santos et al.[[2006| 

t [Gonzalez Delgado & Cid Fernandes 

2010 

) and in dwarf galaxies (e.g., Makarova et al.|2010flGarcfa- 

Benito & Perez-Montero|2012[), it is necessary to improve those 


tests by applying them to more complex systems (i.e. systems 
with complex SFH) with available data of individual stars down 
to the oMSTO in the CMD. 

We have performed such a test using one of the few fully 
resolved bright galaxies where this test can be carried out. 
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the Large Magellanic Cloud (hereafter LMC). In particular, we 
study a region of the LMC bar as a “guinea pig” for which we can 
compare complex SFH’s derived through integrated-light syn¬ 
thesis techniques and through the CMD of the resolved stellar 
population. We emphasize the fact that, in this paper, the SFH 
in the LMC bar region has been derived from the CMD with¬ 
out knowing the results from the integrated spectrum analysis 
and vice versa, making this a “blind test” in order to avoid bias 
toward a common solution by adjusting the fitting parameters. 

Throughout this paper we will consider the SFH as com¬ 
posed by three main functions: one is the star formation rate as a 
function of time, SFR(t); the second is the chemical enrichment 
history (age metallicity relation, AMR); and the third is the ini¬ 
tial mass function (IMF). 

In this paper, we show the LMC bar SFH derived from full 
spectrum fitting techniques applied to its integrated spectrum 
and we compare it to the CMD results from Hubble Space Tele¬ 
scope (HST) data. We have focused our work on the STECKMAP 
results. This code gives more probability to smo oth solutions 


(by means of a penalised fitting algorithm, see |Ocvirk et al. 
|2006b|a| for further information). The smoothness of the recov¬ 
ered SFH makes the comparison with the CMD results easier 
and more sensible. However, we will also compare the results 
obtained using other spectral fitting codes such as ULySS and 
STARLIGHT. In Sect. we explain our target choice as well as 
describe the observations and data reduction procedure to obtain 
the composite CMD and the integrated spectrum. In Sects. |^and 
l^we derive the SFH following both approaches and compare 
the similarities and differences found, with special emphasis on 
how the final solutions are effected by using different STECKMAP 
input parameters. Results coming from other codes (ULySS and 
STARLIGHT) are shown and analyzed in Sect. Some caveats, 
parallel analysis, implications, future work, and main conclu¬ 
sions are discussed in Sects. |6] and |71 


2. Observations and data reduction 

To proceed with the comparison between integrated and resolved 
stellar population approaches in the case of a complex stellar 
population, we selected a field in the LMC bar (see Fig0. The 
centre of the LMC bar is bright enough as to obtain a high quality 
integrated-light spectrum with high S/N, which we did from ob¬ 
servations at the 3.6 m ESO telescope on La Silla using EFOSC2 
( [Alloin et al.|2002| ). It is also sufficiently close and resolved for 
a CMD reaching the oMSTO to be secured with the HST. Such a 
CMD has been already published by |Smecker-Hane et al.| ( |2QQ2] ) 
and various SFHs have been derived by different groups, and 
published in |Skillman & Gallart| ( |2002| ). We present a newly de¬ 
rived SFH with more sophisticated analysis techniques here. 

When the spectroscopic observations were obtained (15 
years ago) the only HST field available with sufficient surface 
brightness to carry out this project was the field observed with 
the WFPC2. Unfortunately, the poor dynamic range of the cam¬ 
era limits our observations of the brightest stars. However, we 
can overcome this limitation with the methodology used to anal¬ 
yse the data (see Sect. 1^. 


2.1. Resolved stellar populations 


The photometry and artificial star tests necessary to recon¬ 
struct the SFH of the LMC bar region were obtained from the 


"HST Local Group Stellar Photometry Archive" (Holtzman et al. 


2006), maintained by J. Holtzmann]^ In particular, we down¬ 
loaded the data of two WFPC2 pointings, namely the u4bll2 
and u4bll5 fields, located at the center of the LMC bar (see 
Fig[T]) and originally observed within the GO program 7382 (PI. 
T. Smecker-Hane). Table details the observing log, including 
the position of the two fields, the integration time in both the 
F555W and F814W filters, and the date of the observations. 

The ^f 555 w~^fuaw) CMD is presented in Fig.[^ 

The photometry reaches down to ~ 6, well below the 

oMSTO. A prominent bright main sequence is visible up to 
^FUAw ~ -1^ indicating that star formation continued until 
very recent epochs in this region. The red giant branch (RGB), 
comprising stars older than ~ 1 Gyr, is also well populated up 
to M/78141V —2. Finally, we also highlight the presence of a 
prominent red clump of centrally He-buming stars {Mp^uw ~ 
-0.5), while the old horizontal branch is barely populated. We 
have marked four regions in the figure that will be discussed in 
Sect.O 


The lack of stars brighter than Mp^uw —1 and Mp^Aw —2, 
along the main sequence and the RGB respectively, is likely due 
to saturation. Indeed, visual inspection of the brightest sources in 
the field confirms that a number of bright stars (~ 200) are satu¬ 
rated in both filters. The region of the CMD affected by satura¬ 
tion will not be used in the SFH derivation. At the metallicity of 
the youngest LMC stars, our saturation magnitude of F814W ~ 
-1 on the main sequence corresponds to the turnoff of a ~ 0.3 
Gyr stellar population. Because the bright massive stars younger 
than this age do not appear in the CMD, this implies that we are 
not using all the possible information to derive the SFH for these 
ages. However, because lower mass stars of the same ages are 
present at fainter magnitudes, the SFH can be still obtained from 
them, under the assumption of a given IMF. 


2.2. Integrated spectrum 


Integrated-light spectra of the LMC bar and a sky field were 
collected at the 3.6 m ESO telescope on La Silla, using EFOSC 
( [Buzzoni et al.|1984| ), on 18-20 December 2000. The J2000 co¬ 
ordinates of the observed target field are: a = 05 : 23 : 17 
and S = -69 : 45 : 42. We used a North-South 5' long and 
1.5” wide slit, which was swept along the East-West direction: 
this allowed us to cover a full area of 2.5' by 5', which approx¬ 
imately matches the WFPC2 pointings (see Figure 1). Four dif¬ 
ferent grisms were used, panning in total the wavelength range 
3500 to 8745A (see Fig. We smoothed the spectra from the 


four grisms to a common final dispersion of 1.9+0.1 A/pix (see 
Fig. [5]). A mean exposure time of 5400 s, divided in three sepa¬ 
rate exposures, was adopted for each grism, both for the target- 
field and also for the sky-field. The sky-field spectrum was taken 
6° North of the target-field position. The sky position was cho¬ 
sen as a compromise between the sky-field being relatively close 
and thus representative of the foreground Milky Way (MW) con¬ 
tamination and sky light present in the LMC bar field, and to 
mimimize the contribution from LMC light. We note, however, 
that with a galactocentric radius of at least 16° ( [Majewski et al. 


|20Q9t |Saha et ar]|2Q10| ), even at 6° from its center, the LMC 
still pres ents a well populated CMD, with m any intermediate- 
age stars ( [Gallart et al.|2008t|Saha et al.|2010| ). 


The 2D-spectra corresponding to the 4 grisms for both 
the target and the sky fields were reduced through standard 


^ http://astronomy.nmsu.edu/holtz/archival/html/lg. 
html 
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RA (J2000) 

5h 22m 30s 23m 00s 



Fig. 1. Left panel: Image of the LMC bar and its surroundings (Credit: John Gleason). The location in the center of the bar of the field studied in 
this paper is indicated. Right panel: Positions of the two WFPC2 fields (red shaded areas), and the 2.5' x 5' area covered by sweeping the slit in 
the East-West direction (blue shaded area), superimposed on a VIMOS B-band image. 


Field 

R.A. 

h m ss.s 

Dec 

deg ' " 

Exp. Time - F555W 
s 

Exp. Time-F814 W 

V 

Date 

u4bll2 

05:22:57 

-69:46:53 

4x500 

2x300-f 2x700 

Nov 27, 1997 

u4bll5 

05:22:55 

-69:42:51 

4x500 

2x300-f 2x700 

Jan 05, 1999 


Table 1. Observing log of the photometric data. 


techniques, using MIDAS and I RAF packages. A full two- 
dimensional wavelength calibration was built up to correct for 
geometrical distortions. The spectra were flux calibrated using 
the spectrophotometric standard EG 21 and LTT 4816; the error 
in flux calibration is estimated to be around 10%. After this step, 
the ID-spectra of the 4 grisms were matched together and the 
quality of the match controlled through the overlapping wave¬ 
length regions. 

To check whether small number statistics in the sampling of 
minority stellar populations in the target fleld (such as the in¬ 
sertion of a few young stars dominating the flnal spectrum but 
with little mass contribution) could lead us to important fluctu¬ 
ations in the SFH, we performed the following test. From the 
integrated-light, spatially resolved 2D-spectra, we extracted two 
series of integrated-light ID-spectra, one for the target-held and 
one for the sky-fleld, with extraction windows of 5' and of 2.5'. 
We ended up with four ID-spectra over the 3500-8745A range 
(FWHM ~ 10 A), two of them corresponding to the integrated 
light in the LMC bar fleld over spatial areas of respectively 
2.5' X 5' (fleld 1) and 2.5' x 2.5' (fleld 2, included in fleld 1) 
and the other two corresponding to the integrated light in the 
sky-fleld with the same extraction windows. At face value, the 
sky-subtracted spectra of the target-held through the two differ- 

^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As¬ 
tronomy (AURA) under cooperative agreement with the National Sci¬ 
ence Foundation. 
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ent extraction windows exhibit very similar characteristics. In 
terms of the shape of their spectra, we scaled the fleld 2 spec¬ 
trum to match the higher flux level of fleld 1 spectrum. The av¬ 
erage difference of both spectra relative to fleld 1 spectrum, over 
the whole spectral range, is 1.6% (i.e. both scaled spectra are al¬ 
most identical). We then derived the SFH of both spectra (fleld 
1 spectr um a nd the scaled, fleld 2 spectrum) using (STECKMAP, 
see Sect. |3.2| ). The residuals between both recovered SFR(t) are 
within our error values (see Sect. |3.2| for further information). 
This analysis indicates that we are properly sampling the stellar 
content in the analyzed fields, in spite of the very small fraction 
of the galaxy’s light present in them. Therefore, we decided to 
analyze in detail the composite stellar population corresponding 
to the spectrum covering the 2.5' x 5' area (40pcx80pc), which 
has a better S/N ratio as a consequence of the larger spatial cov¬ 
erage of an area with intrinsically the same stellar content. The 
corresponding flnal, sky-subtracted spectrum, is shown in Fig.J^ 
The flnal signal-to-noise ratio of this spectrum is ~ 36.3 (per A). 


3. Determination of the star formation history 


We have obtained the SFH using the data described in Sects. |2.1| 
and |2.2[ following well developed and tested methodologies to 
study the stellar conte nt via CMD ( [Aparicio & Gahart| 2004t 
Monehi et al.| 201 0c| ) and spectroscopic analysis ( [Sanchez 


Blazquez et al. ]20TT| |2014t [Seidel et~ar]|2015[ Ruiz-Lara in 


prep.). This study is meant to be a blind test between both ap- 
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Fig. 3. The composite spectrum of the LMC bar field used in the integrated stellar population analysis. Top row: The spectrum is divided according 
to the four different grisms used (different colours). Blue for the grism #7 (3500 - 5185 A), green for grism #9 (4725 - 6730 A), red for grism 
#10 (6275 - 8160 A), and black for grism # 15 (7565 - 8745). For each grism the corresponding sky spectrum, i.e. the spectrum taken in the sky 
field, has been subtracted. Only the LMC bar contribution is shown here. Bottom row: Fully-reduced composite spectrum of the LMC bar (black) 
and fully-reduced, emission-cleaned composite spectrum (red). For the stellar content analysis we use the wavelength range from 3800 to 6800 A 
(see text for details, dashed vertical red lines). The shaded regions of the spectrum are the masked regions (not considered in the fit). 


proaches. Thus, all the results obtained and described in this 
section have been analyzed and studied independently (without 
knowing the results from the other approach) avoiding any pos¬ 
sible bias to a common solution. 


3.1. CMD analysis 


The SFH from the CMD was obtained using the lAC- 


star/MinnIAC/IAC-pop (|Aparicio & Gallart 2004 

Aparicio & 

|Hidalgo||2009[|Hidalgo et al.||2011[|Monelli et al. 

2010cl) suite 


of routines. This SFH derivation is based on the comparison, 
through SiX^ minimization, of the distribution of stars in the ob¬ 
served and in a model CMD. The model CMD was obtained 
from a synthetic CMD computed with I AC-star, after simula¬ 
tion of the observational errors. For the calculation of the syn¬ 
thetic CMD, which contains 9x10^ stars, we assumed a constant 
star formation rate at all ages within 0 and 13.5 Gyr, and fiat 
metallicity distribution between Z = 0.0002 and 0.02. For con¬ 
sistency with the analysis of the integrated spectrum, we used 
the Padova stellar evolution library ( [Girardi et al.|2000bl|Marigo| 


et al.||2008|). Finall y, we assumed a Kroupa initial mass func- 
tion (|Kroupa|2001| ), and a binary fraction of 40% ( [Monelli et al. 


2010c| ). The observational errors were simulated using the infor¬ 


mation from the artificial stars tests. These allow for the inclu¬ 
sion of not only the completeness of the photometry, but also, 
and more importantly, how the photometric errors displace a star 
from its original position in the CMD. The model CMD was di¬ 
vided into simple stellar populations by adopting a basic set of 
age and metallicity bins, whose limits are: age (Gyr) = [0, 0.25, 
0.5, 1.0, 1.5 to 13.5 in steps of 1 Gyr] and Z = [0.0002, 0.0006, 
0.001 to 0.01 in steps of 0.001, 0.015, 0.02]. We adopted a dis¬ 
tance modulus (m-M)o = 18.5 mag and reddening E(B-V) = 0.1 
mag to shift the observed CMD to the absolute magnitude plane. 


Following previous investigations of other dwarf galaxies 


([Hidalgo et al. 

2009 

Monelli et al.|2010c|b| [Hidalgo et al.|2011[ 

[Meschin et al. 

2014 

), the SFH was derived using only the main 


sequence (MS) and sub-giant branch (SGB) regions, avoiding 
the regions where the completeness is less than -50%. Figure 
presents the CMD with the regions adopted for SFH derivation 
{''bundles'') highlighted. The star counts are performed in each 
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^^F555W ^^F814W 


Fig. 2, Mfsi 4 w ^f 555 w ~ ^FSi 4 w CM!D based on ^VFPC2 data. The 
four polygons show the regions used to derive the SFH using the lAC- 
star/MinnlAC/IAC-pop suite of routines (see text for details). 


bundle by dividing it into boxes of different sizes. The finest sam¬ 
pling was used in bundle 1 (0.01 mag and 0.2 mag for the colour 
and magnitude, respectively) and bundle 2 (0.02, 0.1). Bundle 3 
was sampled in (0.2,0.5) boxes and bundle 4 was used as a sin¬ 
gle box. This strategy is an optimal one taking into account that 
i) the physics of main sequence stars is best understood, and ii) 
the largest number of stars in the CMD are located in the lower 
main sequence. Of the 69 209 observed stars, 27 556 are counted 
in the four bundles. 

To take into account possible errors in the distance, redden¬ 
ing, and photometric or model calibration, a number of differ¬ 
ent solutions were derived introducing small shifts in colour 
and magnitude, in a grid of 25 positions within +0.06 mag 
and ±0.15 mag in colour and magnitude, respectively. In or¬ 
der to minimize the effect of the sampling choice, in each po¬ 
sition of the grid 24 solutions were calculated varying the as¬ 
sumed age and metallicity bins of the simple stellar popula¬ 
tions, and moving the boxes inside the bundles. In each position 
of the grid, the 24 solutions are averaged and the correspond¬ 
ing is calculated as the average of the of tho solutions. 

The minimum value of x^ indicates the position in the 

grid where the best solution is obtained, which turned out to be 
A[(M/7555 v^a - M/7814w), A(M/7814w)] = [0.03,0]. 


Figure]^ shows the SFH derived from this analysis. This fig¬ 
ure (blue colours) displays the SFR (t) and Z (t) projections of 
the SFH, as the average of the 24 individual soluti ons calculated 

at the position of the grid where. was obtained (Hidalgo et al. 


2011| ). We should note here that the averaging of the 24 solutions 


has the effect of smoothing the final SFH. Error intervals are cal¬ 
culated from the dispersion of these 24 solutions, plus any solu¬ 
tion in the grid that differs in;^^ by less than Icr. This provides 
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errors equal or in excess to the so-called several solu tions crite¬ 
rion which was shown by |Aparicio & Hidalgo] ( |2009| ) to produce 
reliable estimates of total internal errors. Figure |4| shows that the 
SFR (t) is relatively smooth over the whole time interval, with 
a slight change at around 4.5 Gyr ago an increased SFR(t) at 
later times. This is reminiscent of the period of low star forma¬ 
tion activity followed by a later increase, found by other studies 
in EMC regions located at larger galactocentric distances (e.g., 
[Meschin et al.|20f4^ and references therein). The intermediate- 
age SFR (t) in the EMC bar, however, seems to differ substan¬ 
tially from that of the disc at different galactocentric distances, 
in which the two main periods of star formation, separated by 
an epoch of lower star formation activity, are barely seen in the 
bar and replaced by a much flatter SFR (t). The SFR (t) between 
3.5 and 0.25 Gyr obtained here is rather fiat followed by a de¬ 
crease for stars younger than 0.25 Gyr. The AMR is similar to 
the one reported in previous studies (e.g., [Olszewski et al.|1991t 


Dirsch et al.|2000[|Grocholski et al.|2006] ), showing an exponen¬ 

tial decline from Z ~ 0.0175 at young ages to Z ~ 0.0004 at old 
ages. The metallicity dispersion in the AMR is found to increase 
toward younger ages. 


3.2. Integrated spectrum analysis 


We used the integrated spectrum of the EMC bar (see Sect. |2.2| ) 
to obtain its SFR (t) and AMR employing state-of-the-art tech¬ 
niques taking advantage of the large spectral coverage of our 
data along with some of the most detailed stellar libraries and 
models. 


3.2.1. Emission line cleaning (GANDALF) 

When studying integrated stellar populations it is important to 
take into account any contributions from ionised gas (see Fig. 
S*. To obtain a reliable SFH we must remove such a contribu¬ 
tion, in order to be able to use also those regions affected by 
emission (Hj3, [O iii] A 5007, Ha, etc.) in the stellar population 
analysis. To this end we used GANDALF (Gas AND Absorption 
Fine Fitting, jSarzi et al.|2006|) for the emis sion line removal. In 


this step we made use of thelVazdekis et al.l 

(I20101) models based 

on the MILES library n ([Sanchez-Blazquez 

et al.|2006 

Cenarro 

et al.[[2007[) and computed using the scaled-solar isoc 
Girardi et al.[ ([2000b[). The models are generated fol 

irones of 
lowing a 


Kroupa Universal IMF ( |Kroupa|2001| ). GANDALF simultaneously 
fits the absorption and emission features present in the st udied 
spectrum. It previously runs pPXF (Penalized Pixel Fitting, [Ca^ 
jpellari & Emsellem|2004| ) to obtain the stellar kinematics and a 
combination of spectral templates matching the observed spec¬ 
trum. The stellar kinematics obtained here will be used later on 
(see Sect. [3.2.3] ). Once the absorption spectrum is taken into ac¬ 
count, additional gaussians are included in the fit to obtain the 
kinematics, shapes, and fluxes of the different emission lines. 

3.2.2. Star formation history (STECKMAP) 

The emission-cleaned integrated stellar spectrum was then 
analysed using full spectrum fitting techniques by means of 
STECKMAP (STEllar Content and Kinematics via MAximum 
a Posteriori, [Qevirk et al.[[2006b]a| . The spectrum is fitted us¬ 
ing a Bayesian minimisation method (by means of a penalised 

^ The models are publicly available at http: //miles .iac.es 
^ STECKMAP can be downloaded at http: //astro. u- strasbg. fr/ 
~ocvirk/ 
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to obtain the stellar combination that matches the observed 
spectrum via a maximum a posteriori algorithm. The equation 
to minimise is: 


Qfi = Z, g)) + P^(x, Z, g), (1) 

where v is the modelled spectrum which depends on the age 
distribution (x), the AMR (Z), and the broadening function ac¬ 
counting for the kinematics (g). STECKMAP has the advantage of 
being a non-parametric code, i.e. it does not assume an a pri¬ 
ori shape for the solution (such as stellar age distribution, AMR, 
and the line-of-sight velocity distribution), but gives more prob¬ 
ability to smoother solutions. This smoothness is accomplished 
through the penalisation function (function in equation[^ and 
the function P (a function that gives high values to high oscillat¬ 
ing solutions while small values to smooth solutions of x, Z, or 
g). The penalisation function P^ is dehned as: 


P^(x, Z, g) = PxP(x) + PzP(Z) + yUvP(g). 


( 2 ) 


The parameters Px, Pz, and are the smoothing parameters 
for the Stellar Age Distribution (SAD), the AMR, and the line- 
of-sight velocity distribution (LOSVD), respectively. The differ¬ 
ent values that these parameters can adopt allow us to change 
the smoothness of the accepted solutions. The function P also 
can adopt diff erent shapes (for further information see |Qcvirk| 
et al.||2006b|a| ). This code accounts for the continuum shape of 
the spectrum using a polynomial htting (avoiding sources of er¬ 
ror such as flux calibration or extinction errors). 

STECKMAP outputs allowed us to reconstruct the SFR(t) and 
the AMR that best ht the integrated spectrum (see Fig.|^. Errors 
were computed by means of a series of 25 Monte Carlo simula¬ 
tions. Once STECKMAP has determined the best combination of 
stellar populations to ht the observed spectrum {''best model”), 
we added noise to the "best model” spectrum and ran STECKMAP 
again. This procedure was repeated 25 times. We compute the er¬ 
ror of the stellar mass fraction at each age and metahicity bin as 
the standard deviation of the 25 resulting values. 


3.2.3. Robustness of the STECKMAP results 

There are several input parameters to consider while running 
STECKMAP that might affect the SFH reconstruction (e.g. spectral 
templates, S/N, the smoothing parameters, etc.). In this section 
we check the effect of the various input parameters as well as 
probe the robustness of the STECKMAP results. 

We have explored the entire input parameter space in 24 dif¬ 
ferent tests to choose the be st co mbination of parameters for 
the STECKMAP run (see Sect. |3.2| ). Ah these tests can be inter¬ 
preted as a way of testing the robustness and consistency of the 
STECKMAP results. The key parameters that can affect the hnal re¬ 
sults are: i) the set of spectral templates (models and age range); 
ii) the smoothing parameters; iii) whether the stellar kinematics 
are hxed or htted: although we prefer to hx the stellar kinematics 
in order to stabilise the solution and minimise well reported de¬ 
generacies ( Sanchez-Blazquez et al.||2011| ), we will check how 
the processes of htting simultaneously the stellar content and 
kinematics affects the hnal results. 

Table [^summarizes the main input parameters for each test 
and the corresponding rms. In the tests we explored the whole 
parameter space divided into three main blocks: i) we checked 
the effect of using different stellar models, the age range used 


during the ht, and the simultaneous recovery of the stellar kine 
matics (tests 1 to 8); ii) based on the previous block and our ex 
pertise ( [Sanchez-Blazquez et al.|20ir||2014[ [Seidel et al.|2015 


Ruiz-Lara in prep.) we used the [Vazdekis et al.[ ( [2010[ ) mod- 
els (age-range : 0.063x10^ to 17.8x10^ yr) and hxed the stellar 
kinematics (tests 9 to 16) while exploring the smoothing param¬ 
eters {px and pz) with values ranging from 10“^^ to 10^^; and 
hnahy iii), in the third block (tests 17 to 24) we used the same 
ingredients as in the second block but limited the age range to 
13.5x10^ yr (to match the CMD analysis age range). In ah the 
tests we use a square Laplacian smoothing kernel for the shape 
of the penalisation function, P. 

A visual inspection of the hts, and the quantitative value of 
their rms (~ 0.12 with little dispersion, see Table [^ for ah the 
tests, led us to conclude that there is no easy way to choose 
the combination of input parameters that best suits our data. In 
terms of the reconstructed SFR(t) and AMR shapes, ah of the 
tests show very similar results (compatible within errors), except 
for tests 5, 6, 7, and 8 (see appendix]^ that display the largest 
discrepancies (also in terms of their rms). These tests share the 
property that we ht simultaneously the stellar kinematics, which 
considerably hampers the correct SFH reconstruction. This be¬ 
haviour highlights the overall stability of the STECKMAP solu¬ 
tions, since as the general shape (if not the precise details) of the 
solutions is quite similar among tests. 

As there is no clear set of input parameters to favour over 
the others, we have decided to stick to the following reason¬ 
able set of input parameters: i) We hxed the stellar kinemat¬ 
ics to the values found with pPXF following [Sanchez-Blazquez 


[et al.[ ( [201 1[ ). ii) We matched the age range used in the CMD 
analysis using SSP model templates with ages from 0.063x10^ 
to 13.5x10^ yr. iii) We focused on the results using MILES. In 
Figure we show a comparison between the commonly used 
MILES models ( [Vazdekis et al.[|T010[ ), MILES with the exten 
sion towards youn ger stell ar populations from [Gonzalez D^ 


[gado et aL] ( [20Q5[ ), and the [Bruzual & Charlot[ ( [2QQ3[ ) (BCQ3) 


models. We hnd small differences between models regarding the 
SFR(t). Larger differences are found in the recovered AMR, with 
discrepancies of the order of the errors. In particular, the inclu¬ 
sion of the [Gonzalez Delgado et al.[ ( [2005[ ) models (based on the¬ 
oretical stellar libraries) with ages younger than 63 Myr might 
affect the SFH recovery depending on the real amount of young 
stars present in the observed held. If we do not include those 
models, in principle, the contribution of such stars would be in¬ 
cluded in older bins modifying the recovered SFH shape. Further 
investigation is needed to understand whether the impacts on the 
solutions do (or do not) make sense. 

Bearing this in mind, we show test 24 as an example of 
the recovered SFH. We have used this test just for illustra¬ 
tion purposes, but most of the considered tests could have been 
used without modifying main conclusions. Inspecting Fig.[^(red 
lines) we can see that the LMC bar displays an almost constant 
SFR (t) since its formation until ~ 4 Gyr ago according to the 
analysis of the integrated spectrum. Afterwards, a progressive 
increase in the SFR (t) is found with a peak at ~ 1 Gyr ago fol¬ 
lowed by a drop in the SFR (t) to the present day. The AMR ex¬ 
hibits an exponential increase with old stellar populations show¬ 
ing the lowest metahicity {[MjH] -0.6) and young stars the 

highest metahicities {[MjH] ~ 0.2). We must highlight that this 
exponential behaviour in the AMR shape is mainly caused by 
our choice of smoothing parameters {pz is 10^ for test 24). The 
only difference between tests 17 to 24 is the choice of values of 
the smoothing parameters. As expected (see appendix]^, high 
values of px give smooth SFR(t) shapes (tests 17, 21, and 22) 
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Test 


Models 



Age range (yr) 

Kinematics 

Fx 


Fv 

rms 

1 


Vazdekis et al. 

(2010 

) 

[0.063x10^ 17.8x10^] 

FIX 

10-^ 

102 

_ (**) 

0.1215 

2 

1 

v^azdekis et al. i 

iolo) 

-f 

[0.001xl0^ 17.8x10®] 

FIX 

10-2 

102 

- 

0.1298 


Gonzalez Delgado et al. (2005) 

[0.063x10®, 13.5x10®] 






3 


Vazdekis et al. 

(^0 

) 

FIX 

10-2 

102 

- 

0.1216 

4 

1 

v^azdekis et al. ( 

iolO) 

-f 

[0.001x10®, 13.5x10®] 

FIX 

10-2 

102 

- 

0.1298 


Gonzalez Delgado et al. (2005) 

[0.063x10®, 17.8x10®] 






5 


Vazdekis et al. 

(^0 

) 

FIT 

10-2 

102 

10-2 

0.1670 

6 

1 

v^azdekis et al. ( 

iolO) 

-f 

[0.001x10®, 17.8x10®] 

FIT 

10-2 

102 

10-2 

0.1764 


Gonzalez Delgado et al. (2005) 

[0.063x10®, 13.5x10®] 






7 


Vazdekis et al. 


) 

FIT 

10-2 

102 

10-2 

0.1768 

8 

1 

v^azdekis et al. ( 

ioio) 

-f 

[0.001x10®, 13.5x10®] 

FIT 

10-2 

102 

10-2 

0.1766 


Gonzalez Delgado e 

TaT (2005( 

[0.063x10®, 17.8x10®] 


102 

102 



9 


Vazdekis et al. 

(201(5 

) 

FIX 

- 

0.1218 

10 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 17.8x10®] 

FIX 

10-2 

10-2 

- 

0.1217 

11 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 17.8x10®] 

FIX 

10-® 

105 

- 

0.1224 

12 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 17.8x10®] 

FIX 

10“^ 

10“^ 

- 

0.1217 

13 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 17.8x10®] 

FIX 

10^ 

105 

- 

0.1223 

14 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 17.8x10®] 

FIX 

1015 

1014 

- 

0.1224 

15 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 17.8x10®] 

FIX 

10-15 

10-14 

- 

0.1217 

16 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 17.8x10®] 

FIX 

10-2 

10® 

- 

0.1227 

17 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 13.5x10®] 

FIX 

102 

102 

- 

0.1218 

18 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 13.5x10®] 

FIX 

10-2 

10-2 

- 

0.1218 

19 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 13.5x10®] 

FIX 

10-5 

105 

- 

0.1226 

20 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 13.5x10®] 

FIX 

10-* 

10-* 

- 

0.1218 

21 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 13.5x10®] 

FIX 

10^ 

105 

- 

0.1224 

22 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 13.5x10®] 

FIX 

1015 

1014 

- 

0.1249 

23 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 13.5x10®] 

FIX 

10-15 

10-14 

- 

0.1218 

24 (*) 


Vazdekis et al. 

(2010 

) 

[0.063x10®, 13.5x10®] 

FIX 

10-2 

10® 

- 

0.1225 


Table 2. Set of characteristics and parameters we have used in the different tests to eximane the robust ne ss of the STECKMAP solutions. First column 
is the number of the test. Second column shows the models used in each test, i.e. |Vazdekis et al.| ( [2010| ) or |Vazdekis et al.|p010| )4 )Gonzalez Delgado| 
|et al.|p005| ). Third column represents the age range. In the fourth column we highlight the fact of fitting or fixing the kinematics, stands for 
’smoothing parameter for the stellar age distribution’, pz stands for ’smoothing parameter for the age-Z relation’, and p^ stands for ’smoothing 
parameter for the line-of-sight velocity distribution’. The last column represents the quality of each of the tests by means of its residuals rms 
computed as the mean values of the absolute differences between the data and the fit. (*) Test that we have chosen to compare with the CMD 
results (see text for details). (**) As we are fixing the stellar kinematics, and thus, not fitting it, this parameter is not applicable. 


and high values of pz gives smooth AMR shapes, (tests 17, 19, 
21, 22, and 24) while low values give highly oscillating results. 
In order to easily compare with the rest of solutions we have 
plotted a shaded region corresponding to the average of all the 
solutions using MILES models and ages ranging from 0.063 to 
13.5 Gyr. As can be seen, our chosen test is within that shaded 
region showing again the robustness of the STECKMAP results. 


4. Comparison between CMD and integrated 
spectrum anaiysis 

In Sect. we derived the SFHs in a region of the LMC bar using 
two different techniques. In one technique we used the resolved 
stars in a deep CMD and compared them to theoretical CMDs 
based on theoretical isochrones, in the other we applied full spec¬ 
tral fitting techniques using stellar population synthesis models. 
These analyses were performed independently in a blind test. 
The comparison between the SFH from both approaches will 
then allow us to test how well modern spectral fitting techniques 
recover the characteristics of a complex stellar population. 
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4.1. Star formation rate, SFR(t) 


The upper panel of Fig. shows the SFR (t) in the LMC bar 
region obtained using CMD and STECKMAP typical analyses, as 
discussed in Sect.[^ The overall shape of SFR (t) is qualitatively 
similar in both derivations, with an almost continuous star for¬ 
mation since the earliest epochs till 4 Gyr ago, when small differ¬ 
ences between the SFH obtained from both approaches appear. 
The CMD analysis reveals a rather fiat SFR (t) between ~ 0.25 
and ~ 3.5 Gyr with a clear drop in the star formation during the 
last 0.25 Gyr. The STECKMAP analysis shows a SFR (t) that can 
be described as a gaussian skewed towards older ages with the 
peak around 1 Gyr. Although the youngest populations are bet¬ 
ter sampled in the case of the integrated analysis than in the case 
of the resolved analysis, the decline of SFR (t) at ages younger 
than 0.25 Gyr is consistent in both cases. The lower panel of 
Fig.j^shows the cumulative mass fraction, as a function of time, 
obtained following both approaches. Small differences can be 
noted between the curves representing the mass build-up of this 
region of the galaxy. They are, however, consistent within the 
error bars and hence not significant. 
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Age [Gyr] 


Fig. 4. Comparison between the LMC bar SFH from the CMD and the 
integrated spectrum using STECKMAP (test 24). The three panels show, 
from top to bottom, the normalised SFR(t), the age-metallicity relation, 
and the cumulative mass fraction; with a zoom at young ages for the nor¬ 
malised SFR(t) and the AMR. We plot the envelope of the histograms 
in the main plots in the insets. In order to make a fair comparison a 
normalisation is needed. We normalise the SFR (t) in such a way that 
the sum of the areas of the different rectangles (SFR [M© jyr] x At) is 
1 (upper panel). The AMR plot represents the average metallicity at ev¬ 
ery age bin. Error bars are 1 cr of the resulting distribution of solutions 
from a series of 25 Monte Carlo simulations in the case of STECKMAP 
and 600 different solutions by applying small shifts in the CMD for the 
CMD analysis (see text for details). The shaded regions correspond to 
the mean values and standard deviations of all the solutions using MILES 
and ages ranging from 0.063 to 13.5 Gyr (tests 3, 7, 17, 18, 19, 20, 21, 
22, 23, and 24). 


4.2. Age-metallicity relation 

In Fig. 1^ (middle panel) we show the AMRs from the CMD 
and the STECKMAP study with a zoom at young ages. The over¬ 
all shape and the metallicity range are consistent. Even though 
larger metallicities are systematically found in the STECKMAP 
analysis, they are consistent within the error bars for most of the 
time interval, except the last ^0.25 Gyr. In the spectrum solution, 
there is an upturn from ~ 1.5 Gyr to now, which is not found in 
the CMD analysis, where the metallicity remains constant and 
always below the STECKMAP values. We find the main similari¬ 
ties in the age range between -1.5 and - 3.5 Gyr. At ages older 
than - 3.5 Gyr the AMR from the integrated spectrum displays a 



Fig. 5. Comparison of the different SFHs from STECKMAP, paying spe- 
cial attention t o the choice of models ([Vazdekis et aLl(|2010|),|Vazdekis 
|et al.|(|2010| ) -h [Gonzalez Delgado et al.| ( |2005| ), or |Bruzual & Chariot 
( |2003| )). A zoom at the younger ages for the SFR (t) and the AMR is ap¬ 
plied. The rest of STECKMAP input parameters as in test 24 (see table[5|. 
See caption of Fig. I^for a complete explanation. 


very shallow negative gradient while the AMR from the analysis 
of the CMD shows a steeper negative gradient and lower metal¬ 
licities. A plausible explanation might be found in the choice of 
the smoothing parameters. The best test uses a smoothness pa¬ 
rameter in the AMR of 10^, thus, very smooth solutions for the 
AMR shape, as the one found, are preferred. If we compare this 
solution (Fig. with other solutions using different smoothing 
parameters (see Sect. [3.2. 3] and appendix [B]), we can see that this 
young, metal-rich component disappears for those tests with a 
low value of the AMR smoothness parameter. Thus, the main 
differences regarding the AMR shape (above outlined) might be 
an artifact from the high smoothing parameter imposed on the 
AMR. In fact, tests 17, 18, 20, and 23 (tests with a lower jiz) 
all show quite similar AMR shapes when compared to the CMD 
results. This is also true of the SFR(t) shapes. 


5. Results from other available codes 

For the sake of completeness, we have also used other available 
full spectrum and SED fitting codes (see Sect. [T]). However, as 
this is not meant to be an exhaustive comparison between all 
the available codes for analyzing stellar content from integrated 
spectra, here we analyze the results from those codes we (the 
authors) are more familiar with. The input parameters for these 
tests are as similar to the STECKMAP input parameters as possible. 


5.1. ULySS 

ULyS^ ( Koleva et al. [2009 ) is a full spectrum fittin g code 
which uses Levenberg-Marquardt minimization (from Mark- 


|wardt|2009| ) to fit a linear combination of non-linear parameters. 
It parametrizes the inverse problem as 


^ ULySS can be downloaded at http: //ulyss. univ- lyonl. fr 
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FobsO) = PnO) X {ws VDiVsys, cr, h3, h4) 

( 3 ) 

i=m 

0 ^ Wi Cmpi(ai,a2, /I)), 

i=0 


where is the observed flux at every value of the wave¬ 

length (A); P„ is the multiplicative polynomial and n its degree; 
LOSVD is the stellar line-of-sight-velocity-distribution that de¬ 
pends on kinematics parameters such as the systemic velocity 
(^sys), velocity dispersion (cr), and higher momenta (h3, h4); 
Cmpi are the different components or SSPs dependent on age 
(^i) and metallicity (^ 2 ) with different weights (Wi) computed 
during the fit. 


We have carried out different tests with different input pa¬ 
rameters using ULySS. In a first approach, our model was com¬ 
puted as a sub-set of / = 15 SSPs broadened by a LOSVD from a 
sub-set of the |Vazdekis et al.| ( [2010| ) models (as for the STECKMAP 
analysis). While the ages of the SSPs were fixed between 63 Myr 
and 13.5 Gyr (equally log-spaced) the metallicities were left free 
and they could vary between the limits of the models (-2.3 and 
0.2 dex, the sub-set of models comprises 105 different SSPs). 
We used the emission-cleaned spectrum obtained in Sect. |3.2| and 
fixed the stellar kinematics. We will refer to this test as “SSP” as 
it is based on SSPs as spectral templates. The best fit SFR (t) and 
AMR are shown in Fig. The overall shape of the SFR (t) is 
qualitatively similar to the STECKMAP and CMD reconstructions 
(see Fig.[^, although showing sporadic bursts of star formation, 
which is expected as we are using a combination of SSPs. ULySS 
fails at replicating the AMR at ages younger than ~ 0.5 Gyr; the 
metallicities recovered are lower than those inferred from the 
CMD and the STECKMAP approaches. However, at intermediate 
and old ages (older than 0.25 Gyr) ULySS AMR results are in 
fair agreement with the CMD and STECKMAP ones. 


Although the SSP approach is widely used when recover¬ 
ing the stellar content from integrated spectra, we have also em¬ 
ployed a more complex set of spectral templates, more similar 
to the ones used in the CMD analysis. We made use of the I AC 
MILES webtools|^to create spectral templates via user defined 
SFHs. We decided to use 35 spectra of populations with con¬ 
stant SFR (t) between 0.063-0.178,0.178-0.501,0.501-1.413, 
1.413 - 3.981, and 3.981 - 13.49 Gyr (equally log-spaced in in¬ 
tervals of 0.45 dex) and fixed metallicities ranging from -2.32 
to 0.22 dex [M/H]. These spectral templates may be a better rep¬ 
resentation of the continuous mode of star formation expected 
in real galaxies than discrete SSPs. We will name these new 
spectral templates “complex SPs” (complex stellar populations). 
Thus, we use these i = 35 “complex SPs” (Cmpt), following 
the same approach outlined above (input spectrum, wavelength 
range, fixed stellar kinematics, etc). We will refer to this test 
as “constant SFR”. The results are shown in Fig. [ 7 ] The use of 
complex SPs considerably improves the recovery of the shape of 
the SFR (t). The main discrepancies are found at ages younger 
than 0.5 Gyr, with an absence of populations with ages between 
0.178 - 0.501 Gyr and an excess contribution in the youngest bin 
(0.063 - 0.178 Gyr). The AMR is similar to the one we obtain 
using SSPs, with the youngest ages showing too low metallici¬ 
ties. 



Age [Gyr] 


Fig. 6. Top panel: SFH using ULySS with Vazdekis et al. ( |2010| ) models 
younger than 13.5 Gyr (“SSP” test). Bottom panel: Age-metallicity re¬ 
lation from the ULySS results. SFH: Note that the representation of the 
SFH is different that in the case of STECKMAP (see Fig.[^. In this case 
we are using SSPs, with no smoothing applied, and thus, vertical lines 
are chosen to show the mass contribution of the different SSPs instead of 
a bar plot. However, for a fair comparison, a similar normalization has 
been applied. Horizontal lines represent the At used in this case (com¬ 
puted based on stellar population models used). For further information 
see Fig. AMR: Green points represent the age and metallicity of the 
different single stellar populations with a non-zero weight in the fit. The 
point size is proportional to such weight. 


5.2. STARLIGHT 


2005| > FI 


We have also used STARLIGHT ( |Cid Fernandes et al 
to obtain the SFH from the observed integrated spectrum. 
STARLIGHT is a SED fitting code that tries to match the observed 
spectrum (continuum and spectral features) by means of a com¬ 
bination of stellar models. It also fits the reddening and stellar 
kinematics simultaneously. 


The details of how STARLIGHT works are given in |Cid FeF] 
[nandes etar| ( |2004| ) and jCid Fernandes et al.| ( |2005'] ). STARLIGHT 
mainly requires as input an observed spectrum, a configuration 
file, a mask file, a set of base spectra (templates), and a red¬ 
dening law. Essentially, the code then tries to obtain the SFH, 
reddening, and stellar kinematics by the minimisation of a;^^: 



( 4 ) 


^ http://miles.iac.es/pages/webtools/ 
get-spectra-for-a-sfh.php 
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^ STARLIGHT can be downloaded at http://astro.ufsc.br/ 
starlight/ 
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Fig. 7. Same as Fig.[^ but using as spectral templates a set of 35 spec¬ 
tra generated from a constant SFR (t) and single metallicity (see text 
for details). However, as in this case we are not using SSPs as spectral 
templates, we use again a similar representation to the one used with 
STECKMAP (see Fig.Q. 


where is the observed flux, is the modelled flux, and 
is the weight (0 for masked regions) at the wavelength A. 
STARLIGHT normalises the observed spectrum and the spectral 
templates at a given wavelength, in our case we use the window 
between 5590 and 5680 A for the observed spectrum, and the 
flux at 5635 A for the template normalization. The expression 
that this code uses for the modelled spectrum is as follows: 




Mx = Xjbj^xrx] ® G(v*, cr*), 


(5) 


i=i 


where M^o is the synthetic flux at the normalization wave¬ 
length; is the modelled flux at A\ xj is the weight of the ele¬ 
ment of the set of base spectra; is the normalised reddened- 

spectrum for such component where 

extinction term. This first part of the expression accounting for 
the stellar content and the reddening is convolved ((8)) by the 
LOSVD of the stellar component (G(v^,(T^)). We made use of 
the LMC reddening law provided by the STARLIGHT package 
and presented in [Gordon et al.| ( [2003| . 


We have analysed the LMC bar integrated spec trum fol low- 
ing a recipe similar to the one discussed in Sects. [T2| and [5.1 [ 
this test will be called “SSP” as it will be based o n SSPs as spec¬ 
tral te mplates. We used the same set of models p|( |Vazdekis et al.j 
2010| up to 13.5 Gyr) as in previous cases, removed the emis¬ 
sion line contribution using GANDALF, fixed the stellar kinemat¬ 
ics to the pPXF values, and run STARLIGHT to match every spec¬ 
tral feature in the wavelength range from 3800 to 6800 A. We 
masked transition regions between grisms and sky features, but 
not the emission lines as they have meaningful information af¬ 
ter applying GANDALF. We And that the results from STARLIGHT 
show some important differences when compared to the CMD 
and STECKMAP results (see Fig.[^. The recovered SFR (t) shows 
some episodic bumps of star formation younger than 2 Gyr and 
a predominant old stellar population (older than 10 Gyr), while 
an intermediate population between 2 and 10 Gyr is not found. 
The main discrepancies regarding the AMR are found at ages 
younger than 0.5 Gyr. As in the case of ULySS, the discontinu¬ 
ity in the recovered SFR (t) is a direct consequence of the use of 
SSPs as spectral templates. 

These discrepancies encouraged us to carry out a set of 36 
tests modifying different input parameters and procedures to 
check the reliability of the results, i) We tested the effect of 
changing the set of model templates and the age range, ii) To re¬ 
duce the degrees of freedom in the fit we also inspected the effect 
of imposing an a priori AMR from jCarrera et al.| ( |2008| ), using 
a carefully selected set of spectral templates from the different 
models, iii) Considering that STARLIGHT also fits the continuum 
shape of the spectrum, we have also tried to obtain stellar con¬ 
tent with this code avoiding the emission line removal step with 
GANDALF, masking the emission lines instead in order to test if 
the discrepancies are caused by a bad emission line subtraction, 
iv) We also allowed STARLIGHT to fix or fit the stellar kinemat¬ 
ics. 


We obtained a wide variety of solutions from these different 
tests. A clear bi-modality is found in the SFR (t) shape using the 
Vazdekis et al.| (|201Q| ) and jVazdekis et al.| ( |2010| ) plus [Gonzalez 


Delgado et al.[ ( 2005[ ) models, with contributions from young and 


old populations (younger than 2 Gyr and older than 8 Gyr) and 
a lack of an intermediate component. However, the use of BC03 
models results in a more spread-out, or smoother SFR (t) with 
contributions at all ages. Imposing an observed AMR strongly 
restricts the number of templates used, and gives poor fits. Fixing 
the kinematics gives slightly better results than fitting for the 
kinematics. Using GANDALF or masking the emission lines plays 
a minor role in the recovery of the stellar content. 

In addition, as we did with ULySS (see Sect. [5.1[ ), we per¬ 
formed a further test where we used complex SPs, rather than 
SSPs, as spectral templates (see Sect. 1^. We name this test 
“constant SFR” as previously done with ULySS. The results of 
this test are shown in Fig. Using complex SPs we have been 
able to obtain a closer approximation to the CMD results re¬ 
garding the SFR (t) shape. However, we And a lack of stars with 
ages ranging from 0.501 - 1.413 Gyr and an excess of stars in 
the youngest age bin (0.063 - 0.178 Gyr). The obtained AMR 


^ In the case of STARLIGHT, as well as for STECKMAP, we use the entire 
set of models, [ages (Gyr)]x[M/H] = [0.0631, 0.0708, 0.0794, 0.0891, 
0.1000, 0.1122, 0.1259, 0.1413, 0.1585, 0.1778, 0.1995, 0.2239, 

0.2512, 0.2818, 0.3162, 0.3548, 0.3981, 0.4467, 0.5012, 0.5623, 

0.6310, 0.7079, 0.7943, 0.8913, 1.0000, 1.1220, 1.2589, 1.4125, 

1.5849, 1.7783, 1.9953, 2.2387, 2.5119, 2.8184, 3.1623, 3.5481, 

3.9811, 4.4668, 5.0119, 5.6234, 6.3096, 7.0795, 7.9433, 8.9125, 

10.0000, 11.2202, 12.5893, 14.12]x[-2.32, -1.71, -1.31, -0.71, -0.4, 0.0, 
0 . 22 ] 
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Fig. 8. Top panel: SFH from STARLIGHT with Vazdekis et al. ( 2010| ) 


models younger than 13.5 Gyr. Bottom panel: Age-metallicity relation 
from the STARLIGHT results. SFH: Note that the representation of the 
SFH is different that in the case of STECKMAP (see Fig.[^. In this case 
we are using SSPs, with no smoothing applied, and thus, vertical lines 
are chosen to show the mass contribution of the different SSPs instead of 
a bar plot. However, for a fair comparison, a similar normalization has 
been applied. Horizontal lines represent the At used in this case (com¬ 
puted based on the base of stellar models used). For further information 
see Fig.|^ AMR: Green points represent the Age and Metallicity of the 
different single stellar populations with a non-zero weight in the fit. The 
point size is proportional to their weights. 
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Fig. 9. Same as Fig.[^ but using as spectral templates a set of 35 spec¬ 
tra generated from a constant SFR (t) and single metallicity (see text 
for details). However, as in this case we are not using SSPs as spectral 
templates, we use again a similar representation to the one used with 
STECKMAP (see Fig.Q. 


is very similar to the one recovered using SSPs as spectral tem¬ 
plates (again, similar to the ULySS results), displaying the same 
issues at the youngest ages. 


6. Discussion 

The LMC bar is an appropriate astronomical object to perform 
a comparison between the SFH obtained for a composite stellar 
population by means of the information from its resolved stars 
and from its integrated light. The reason is twofold: i) we can 
obtain a CMD reaching the oMSTO using HST data, and ii) its 
high surface brightness allows us to observe a high quality in¬ 
tegrated spectrum. By analyzing the CMD and integrated spec¬ 
trum of the LMC bar we derived and compared its SFHs. In this 
work we mainly focused on the comparison between CMD anal¬ 
ysis and full spectrum fitting analysis using STECKMAP ( |Qcvirk| 
|et al.|2006b|a| ). In principle, the stellar populations of the LMC 
bar represent a challenging case for inversion codes analyzing 
integrated stellar spectra because of its complex nature and rich 
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stellar diversity. However, the agreement in the SFR(t) and the 
AMR between the CMD and the integrated spectrum analysis 
strongly supports the use of state-of-the-art full spectrum fitting 
codes in deriving the SFH of complex stellar systems. The re¬ 
sults of our “blind tests” are that the results of both approaches 
are consistent when trying to recover the SFH of the same ob¬ 
ject. However, during our analysis we identified several issues 
that should be further studied. 

The most striking difference between the STECKMAP and 
the CMD results appears in the derived AMRs. Both AMRs 
show a metallicity monotonically increasing with time, but the 
AMR derived from the CMD starts at lower metallicities (^ - 
1.2 dex) and does not reach solar metallicity at young ages. The 
AMR obtained with STECKMAP starts at [M/H] ^ -0.75 dex at 
old ages, remains quite fiat for most of the time range and in¬ 
creases steeply in the last ^ 2 Gyr, reaching metallicity slightly 
over solar at the present time. The overall SER(t) shape also 
displays some differences for populations younger than 4.0 Gyr. 
The SFR (t) is almost flat between ~ 0.25 and ~ 3.5 Gyr from the 
CMD analysis, whereas from the STECKMAP analysis we see a 
skewed gaussian towards older ages with the peak at ~ 1 Gyr. 
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There are several possible explanations for these differences in 
SFR (t) obtained from the two methods: 

a) The CMD analysis produces precise results when a deep 
CMD reaching to below the oMSTO is used (as in this case). The 
methodology described in Sect. |3.1| has been extensively tested 
by recovering SFHs using mock CMDs with known stellar con- 
tent (e.g. |Aparicio & Hidalgo [2009 [ [Monelli et al.|201Qa|ct Hi-| 


jdalgo et al.||2011| )r However, these analyses are not error-free 

Differences in the solutions are found due to slight changes in 
the analysis method (e.g. way of dealing with errors, bundle def¬ 
inition, minimisation algorithm). Also, tests with mock stellar 
populations reveal that an age-dependent smoothing of the SFH 
features is produced, with age resolution worsening at old ages, 
due to the limitations in the age resolution intrinsic to the method 
(e.g. see hgure 8 in [Hidalgo et al.|2QTT] ). Having a large number 
of stars in the CMD is important for a reliable solution, and we 
cannot exclude the possibility that the shape of the SFR(t) at 
young ages is not well recovered due to small number statistics. 
This may be aggravated by the fact that the brightest, youngest 
stars in the main sequence are saturated, and therefore do not 
appear in the CMD. Even though the SFR (t) at the correspond¬ 
ing ages may still be recovered from lower mass stars, which are 
fainter on the main sequence, some information is obviously lost. 
Thus, further investigation using a CMD taken over a larger area, 
and covering the whole magnitude range would be of interest. 

b) There could be a lack of young stars in the held where the 
integrated spectrum was observed. This might produce a dehcit 
of hux coming from the youngest populations, and therefore 
might result in the failure of STECKMAP to reproduce the AMR at 
such young ages. However, we note that STECKMAP, STARLIGHT 
and ULySS ah hnd some contribution at very young ages (see 
Fig.ig. In addition, a visual inspection of the observed spectrum 
shows some helium absorption lines (He I/I3819 and He I/I4922), 
which are signatures of the presence of young, hot stars. There¬ 
fore, it does not seem that a dehcit of signal from young stars 
is the cause of the differences in the metallicity inferred for the 
very young population. 

c) The AMR shape is especially affected by the smoothing 
penalty function in the AMR (/i^). In particular, the shape of our 
recovered AMR using STECKMAP (see Fig. is very smoothed 
as a consequence of imposing a large smoothing parameter (//^ 
= 10^). The discrepancies between this AMR and that from the 
CMD analysis in the younger and older edges of our age range 
could be due to the effect of this parameter. In fact, tests with 

above 10^ currently present supersolar stars and primordial 
metahicities greater than [MjH] = -1.0. Other tests with different 
smoothing parameters for the AMR (see appendix]^ are better 
able to reproduce the AMR shape derived in the CMD analysis 
while still providing a good ht for the SFR(t). 

d) As shown in Sect. |3.2.3[ different input parameters give 
similar SFH results and very accurate hts (see Fig. [T^ and ta- 
ble[^ when STECKMAP is applied to a high quality spectrum. The 
overall shape of the recovered SFH is reasonably consistent for 
the different tests although slight differences do arise between 
them. Therefore, we conclude that the shape of the recovered 
SFH from high quality spectra is well reproduced regardless of 
the input parameters. 

Despite these small discrepancies, the similarities between 
the CMD and the STECKMAP results are very reassuring. The 
manner in which STECKMAP deals with the intrinsically ill- 
posed inversion problem (regularization through a penalized;^^) 
has proven to be very powerful. The smoothed solutions from 
STECKMAP are perhaps physically “sensible”; star formation in 
galaxies is a complex mechanism, and a consequence of continu¬ 


ous galaxy evolution. In STARLIGHT and ULySS no such smooth¬ 
ing is implemented, and thus, the use of SSPs as spectral tem¬ 
plates naturally leads to a SFR (t) characterized by discrete bursts 
of star formation. Although the use of SSPs as a base for recon¬ 
structing the SFH of galaxies is a good hrst approximation, these 
SSPs are idealized realizations of bursts of star formation that are 
not necessarily expected in nature. The inclusion of some kind 
of smoothing techniques in these codes might eliminate these 
discontinuities. 

Since the exact determination of the shape of the SFR (t) is 
difficult in STARLIGHT and ULySS if SSPs are employed, we 
used a set of spectral templates computed assuming populations 
with a range of ages (complex SPs). Although the shape of the 
AMR is very similar irrespective of whether we use SSPs or 
complex SPs, the shape of the SFR (t) obtained using this second 
approach is closer to the CMD and STECKMAP SFR(t) shapes. 
The use of complex SPs seems to in some sense mimic the effect 
of imposing a smoothing on the solution (see Figs. an d|3 and 
it is preferred to the use of SSPs as it allows us to obtain some 
intermediate contribution that otherwise is impossible to detect. 
This suggests that these kind of spectral templates may be better 
in order to study the stellar content of external complex systems 
(e.g. galaxies, not star clusters). 

We have also evaluated whether, in spite of not always pro¬ 
ducing consistent SFR(t) shapes, the different methods are able 
to reproduce the relative contribution in mass of young, inter¬ 
mediate, and old components (the age intervals for this test are 
dehned to match the analysis a la Bica, see Appendix |A|). Table 
[3] shows the contributions to the total stellar mass for different 
age ranges, as inferred using ah the analysed approaches. The 
outcomes from the CMD, STECKMAP, STARLIGHT and ULySS 
(“SSP” and “complex SSP” tests) display roughly similar per¬ 
centages of stars in each age range, although slightly larger dis¬ 
crepancies are found in the STARLIGHT analysis. 

A common issue in the ULySS and STARLIGHT results is the 
excess of young stars (see Figs. muis and[^. As discussed 
in Sect. |2.1| some bright stars (F814W ~ -1) are saturated, and 
thus, not taken into account in the SFH recovery from the CMD. 
This saturation magnitude roughly corresponds to the main se¬ 
quence turnoff of a ~ 0.3 Gyr old population. However, we con¬ 
sider in our calculation other stars with the same age that are 
not saturated (lower mass and therefore, fainter stars), so we can 
still recover the information at the youngest ages. As this satu¬ 
ration affects the brightest (and youngest) stars which might be 
dominating in the light integrated spectrum, the excess of young 
stars found with ULySS and STARLIGHT may be a consequence 
of this limitation. We note that this issue is somehow observed 
in the STECKMAP results with low smoothing parameters. This 
should be further investigated with the comparison with other 
local group galaxy using similar analysis to this work. 

In Fig. [T^ we analyze each spectral ht (STECKMAP, ULySS 
and STARLIGHT) in order to assess where the observed discrep¬ 
ancies discussed above (SFH and age percentages) may come 
from. At hrst glance the three codes seem to properly ht the 
observed spectrum. The rms of the residuals is 0.12 for the 
STECKMAP ht, 0.21 for the STARLIGHT ht and 0.18 for the ULySS 
ht (both “SSP” and “constant SFR” tests). A careful inspection 
shows that STECKMAP and ULySS are able to better reproduce 
some individual spectral features than STARLIGHT as a conse¬ 
quence of the polynomial htting. Figure [m shows a more de¬ 
tailed comparison between the STECKMAP and STARLIGHT best 
models (“constant SFR”), i.e. the spectrum corresponding to the 
recovered stellar content. We note some important differences in 
the shape of the continuum (wave-like features in the residual 
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Approach 

Young 

(%) 

Intermediate 

(%) 

Old 

% 

CMD 

4.3 

49.0 

46.7 

STECKMAP 

4.2 

49.4 

46.4 

ULySS“SSP” 

3.9 

51.4 

44.7 

ULySS “const. SFR” 

2.4 

51.3 

46.3 

STARLIGHT “SSP” 

9.0 

54.9 

36.1 

STARLIGHT “const. SFR” 

7.7 

33.1 

59.2 

Bica (*) 

2.0 

10.0 

88.0 


Table 3. Percentages of the total stellar mass in three different stel¬ 
lar sub-populations: young (younger than 0.5 Gyr), intermediate (older 
than 0.5 Gyr and younger than 5.0 Gyr), and old (older than 5.0 Gyr). 
In the cases of ULySS and STARLIGHT we show the results from the 
SSP and the mock spectral template approaches. We include the Bica 
(*) analyses for historical comparison. For information about the Bica 
analysis see Appendix [A| 



a(A) 


Fig. 11. Comparison between the STECKMAP and STARLIGHT (using 
spectral templates mimicking constant SFR) best models. Top panel: 
The STECKMAP best model is plotted in red. The STARLIGHT best 
model is plotted in blue. Bottom panel: Residuals of both best models 
(STECKMAP - STARLIGHT). The shaded regions of the spectrum are the 
masked regions (not considered in the fits). Vertical, dashed red lines 
are located at the wavelengths of the main stellar absorption features. 
The residuals are in units of 10“^^ erg“^ cm“^ 


plot) as well as differences in some specific absorption features. 
Among the features with larger discrepancies we can highlight 
He, Hp, the Can doublet (3933, 3969 A), the magnesium fea¬ 
ture at 5175 A, and the sodium absorption line (5892 A). Every 
stellar feature is deeper in the STECKMAP best model than in the 
case of the STARLIGHT fit except for the sodium feature. These 
differences affect considerably the recovered stellar content. 

Based on the work presented here we find that, although the 
techniques for the analysis of unresolved stellar populations give 
good results when compared to the CMD reconstruction, there is 
still room for improvement. For example, a wider range of spec¬ 
tral coverage might help to overcome the issues at young ages 
and to better constrain the recovered SFH. Few empirical stel¬ 
lar libraries are available in the ultraviolet (150 to 333 nm) and 
near-infrared (up to 2 microns) and thus, theoretical SSPs limited 
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to the optical range are used (e.g., 

Riffel et al. 2009 Storchi- 

|Bergmann et al.||2012 

Martins et a 

L||2010 2013). In addition, 


we find that a successful way of overcoming the ill-conditioned 
problem of the reconstruction of the stellar content from an ob¬ 
served spectrum is by imposing smoothed solutions (which have 
more physical meaning), or by using spectral templates com¬ 
posed of complex SPs. We suggest that this could be a basis for 
improvement in the next generation stellar population synthesis 
codes from integrated spectra. 


7. Conclusions 


In this paper we have studied the SFH and the AMR of the 
EMC bar using two distinct approaches. We compared the results 
from the reconstruction of an observed CMD reaching the oM- 
STO, obtained from WFPC2@HST data, with the results from 
the spectral fitting of an integrated spectrum taken at the 3.6 m 
ESO telescope on Fa Silla, using EFOSC2. Due to the relatively 
high surface brightness of the EMC bar, and its close proxim¬ 
ity, the chosen field is an unique one in which both an accurate, 
deep CMD and a high quality spectrum can be obtained. We have 
applied state-of-the-art models, stellar libraries, isochrones, and 
codes following each approach in order to derive the SFR (t) and 
AMR of this bar field. Different codes (STECKMAP, ULySS, and 
STARLIGHT) have been used to recover the stellar content from 
the integrated spectrum and the results compared to the CMD 
analysis. The analysis of the integrated spectrum using each code 
has been performed in a consistent manner, and in an effort to 
avoid any biases in the solutions, the CMD analysis was per¬ 
formed independently from the spectral analysis. 

The best agreement between the integrated spectrum analysis 
and the CMD analysis was found using STECKMAP, the only full 
spectrum fitting code that we tested with a penalization imple¬ 
mented. STECKMAP produces SFR (t) and AMR in good agree¬ 
ment with those obtained from the CMD. All the spectral fitting 
codes used in this study are normally used with SSPs as spectral 
templates. ULySS and STARLIGHT do not use any penalization 
and, as a consequence, solutions dominated by episodic bursts 
are derived if SSPs are used as spectral templates. ULySS is able 
to reproduce the overall shape of the SFR(t) but with ’bursts’ 
of star formation. However, this code fails at reproducing the 
AMR, especially at young ages. STARLIGHT is able to approx¬ 
imately reproduce the AMR (except at the youngest ages), but 
has problems with the shape of SFR (t), especially at intermedi¬ 
ate ages where no equivalent contribution in the CMD is found. 
We have been able to improve these results (both with ULySS and 
STARLIGHT) by using a set of complex spectral templates con¬ 
structed adopting a constant SFR (t) in bins of log-spaced ages 
instead of SSPs. This suggests that “complex SPs”, rather than 
simple SP spectral templates should be preferred when analyzing 
the stellar content of composite stellar systems. 

This is the first time that the results of these two different ap¬ 
proaches for studying stellar populations in galaxies have been 
compared for an object with a complex SFH, and for which 
a CMD reaching the oMSTO could be obtained. Such stud¬ 
ies are of crucial importance in order to test recent advances 
in the field, especially in the analysis of the integrated stellar 
populations. Important ongoing and upcoming projects such as 
CAFIFA dSanchez et aTl[MT^ , SAMI ( |Croom et aI.1[20T2l ) or 
MANGA will make use of these techniques to study the forma¬ 
tion and evolution of galaxies. 

In future work we plan to expand our analysis to a sample 
of Focal Group dwarf galaxies bracketing a range of properties. 
This will allow us to identify the impact that different factors 
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Fig. 10. Top panels: Comparison between the observed spectrum and the fit from STECKMAP (left), STARLIGHT (middle), and ULySS (right). For 
STARLIGHT and ULySS we are using the “constant SFR” approach. Black: Fully-reduced, emission-cleaned composite spectrum of the LMC bar. 
Red: Fit from the different codes. Bottom panels: Residuals of those fits computed as observed - best model. For the stellar content analysis we 
use the wavelength range from 3800 to 6800 A (see text for details, dashed vertical red lines). The shaded regions of the spectrum are the masked 
regions (not considered in the fits). The residuals are in units of 10“^^ erg"^ cm"^ 


(such as different SFHs, fractions of young vs. old stars, the ex¬ 
istence of a blue horizontal branch or blue straggler stars) may 
have on the results. Not only will this comparison identify where 
full spectrum fitting techniques may fail and where they suc¬ 
ceed, but it will also inform us on the spectral ranges and fea¬ 
tures which are more likely to improve these techniques. Under¬ 
standing the limits of the reliability of SFHs obtained from inte¬ 
grated spectra provides an important basis for the understanding 
of galaxies at low and high redshift. 
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Appendix A: Analysis using a base of integrated 
cluster spectra (Bica method) 

To historically link previous methods to study the stellar con¬ 
tent in galaxies with modern SED and full spectrum fitting tech¬ 
niques, we analyze the data via a spectral population synthesis 
technique originally developed by |Bica|fl988 ), later updated by 
[Schmitt et al. \\996) . This method aims at reproducing the ob¬ 
served equivalent widths (Ws) and the continuum ratios (Cs) 
using the integrated-light spectra of an ensemble of star clus- 
ters with different ages and metallicit ies (e.g., |Bica & Alloin 


|1986a|bl \TWf\ IBica et al.||T988| |T^ . The Ws and Cs values 
from the cluster base are built in a grid parametrized by the age 
and the metallicity, and extrapolated in the case of high metal¬ 
licities unreachable through observations. In the current analysis, 
we have used eight components to map the age-metallicity plane, 
with a constraint on their metallicities to be solar and subsolar. 
Indeed, in a low mass galaxy such as the EMC, we do not expect 
to find stellar components with metallicity above solar ( [Pagel &| 
Tautvaisien^[1998[ ). The base elements used in the a nalysis are 


listed in table [A. 1[ The results are shown in Eig. [A.l| 

Erom the population analysis above it is clear that the 1 Gyr 
and 100 Myr components are very important in flux (see Eig. 
[A.l[ top panel). It is also important to compute how much they 
represent in terms of mass fractions. Eor such purposes we em¬ 
ployed a flux-mass transformation method ( [Bica et ar][1988| ). 
This method uses different mass to V-light ratios (M/Ly) related 
to each age component. It also takes into account metallicity 
effects among old star clusters. We show in Eig. [A.l[ (bottom 
panel) the mass distributions respectively for each component. 
We can see that the 1 Gyr and 100 Myr components are less than 
10 % and 5 % in mass respectively. 

We can interpret those results as star formation along three 
different age ranges (less than 0.5 Gyr, 0.5 - 5 Gyr, and 5 - 
14 Gyr, see Table [A.2| ). If we consider that the mass at ev¬ 
ery age range has been formed uniformly and we compute a 
“pseudo-SER” dividing the mass alive percentage by the age 
width for each population, that leads us to recent star forma¬ 
tion of 4.4 Gyr“\ followed by a drop (2.2 Gyr“^), with a higher 
“pseudo-SER” at older ages. These results are away from the 
STECKMAP and CMD results (see table 3), essentially because 
of a lack of time resolution at old ages (see table |^. 
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Range 

^alive 

(Gyr) 

(%) 

5-14 

88 

0.5-5 

10 

<0.5 

2 


Table A.2. Results for the Bica 


analysis. MaUve 


percentage in mass that is still present as stars. 



Fig. A.l. Synthesis results using the Bica analysis. Top panel, flux frac¬ 
tion of the different base elements. Bottom panel, mass alive fraction of 
the different base elements. 


3 Myr 10 Myr 100 Myr 1 Gyr 

10 Gyr 

[M/H] 

8 7 6 4 

1 

0.0 

5 

2 

-1.0 


3 

-2.0 


Table A.l. Base elements in the age x metallicity plane for the |Bica 
( |1988| ) analysis. E(B-V)/ = 0.0. The different base elements are identi¬ 
fied by numbers running from 1 to 8. 


Appendix B: Recovered Star Formation Histories 
for the STECKMAP tests 

We include in this ap pendix the recovered SFH in the tests de¬ 
scribed in Sect. l3.2.3l We show the stellar content for the 24 tests 
in a similar as in Fig.|^ As already outlined in the main body of 
the paper, we must highlight the consistency between different 
tests except in some extreme cases. 


represents the 
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Fig. B.l. Comparison between the LMC bar SFH from the CMD and the integrated spectrum using STECKMAP. Left: Test 1; Right: Test 2. For 
further information see Fig.|^ 




Fig. B.2. Comparison between the LMC bar SFH from the CMD and the integrated spectrum using STECKMAP. Left: Test 3; Right: Test 4. For 
further information see Fig.|^ 
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Fig. B.3. Comparison between the LMC bar SFH from the CMD and the integrated spectrum using STECKMAP. Left: Test 5; Right: Test 6. For 
further information see Fig.|^ 




Fig. B.4. Comparison between the LMC bar SFH from the CMD and the integrated spectrum using STECKMAP. Left: Test 7; Right: Test 8. For 
further information see Fig.|^ 
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Fig. B.5. Comparison between the LMC bar SFH from the CMD and the integrated spectrum using STECKMAP. Left: Test 9; Right: Test 10. For 
further information see Fig.|^ 




Fig. B.6. Comparison between the LMC bar SFF[ from the CMD and the integrated spectrum using STECKMAP. Left: Test 11; Right: Test 12. For 
further information see Fig.|^ 
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Fig. B.7. Comparison between the LMC bar SFH from the CMD and the integrated spectrum using STECKMAP. Left: Test 13; Right: Test 14. For 
further information see Fig.|^ 




Fig. B.8. Comparison between the LMC bar SFF[ from the CMD and the integrated spectrum using STECKMAP. Left: Test 15; Right: Test 16. For 
further information see Fig.|^ 
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Fig. B.9. Comparison between the LMC bar SFH from the CMD and the integrated spectrum using STECKMAP. Left: Test 17; Right: Test 18. For 
further information see Fig.|^ 




Fig. B.IO. Comparison between the LMC bar SFH from the CMD and the integrated spectrum using STECKMAP. Left: Test 19; Right: Test 20. For 
further information see Fig.|^ 
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Fig. B.ll. Comparison between the LMC bar SFH from the CMD and the integrated spectrum using STECKMAP. Left: Test 21; Right: Test 22. For 
further information see Fig.|^ 




Fig. B.12. Comparison between the LMC bar SFH from the CMD and the integrated spectrum using STECKMAP. Left: Test 23; Right: Test 24. For 
further information see Fig.|^ 
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